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Notation 


CD  :  Drag  Coefficient 

CT  s  Lift  differential  coefficient 

L(x 

C  :  Static  moment  differential  coefficient 

in 

a 

C£  :  Resultant  force  of  perturbing  influences 
d  :  Diameter 

F  :  Force 

g  :  Acceleration  of  gravity 

I  i  Moment  of  inertia  around  lateral  axis  of  projectile 

Ix  :  Moment  of  inertia  around  longitudinal  axis  of  projectile 

nt  :  Mass 

M  :  Moment 

Me  :  Resultant  moment  of  perturbing  influences 
p  :  Angular  velocity  around  Xg  axis 

q  :  Angular  velocity  around  yfi  axis 

r  :  Angular  velocity  around  Zg  axis 

rt  :  Distance  between  center  of  gravity  and  nozzle  throat 

v  ;  Ve loci ty 

w  :  Wind  velocity 

x,y,z  :  Ground  axis  coordinates 

a  :  Angle  of  attack 

6  :  Side  slip  angle 

0,5  :  Angles  defined  in  Chapt.  3  (2) 

♦  ,£  t  Angles  defined  in  Chapt..  3.  (3) 

n  t  Angle  rotated  around  Xg  axis 


y  Factor  angle  of  trim  of  canard 

v  :  Trim  angla  of  canard 

ir  :  Factor  angla  caused  by  attitude  of  projectile 
p  :  Density  of  air 

Suffix  : 

a  :  Wind  axis  coordinate 

B  :  Body  axis  coordinate 

v  :  Velocity  axis  coordinate 

x,y,s  :  Rectangular  axis 
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1.  IMTBOPPCTIOH 


Solid  body  equations  to  calculate  tha  trajectory  of  a  projactlla 
which  la  flying  with  thrust  of  ramjet  under  hypersonic  Mach  ntaber  are 
forwulatad*  The  projectila  is  gun  launched.  This  projectile  has  an  axl- 
■ywnatrlc  configuration  and  is  fin  stabilised  with  four  pop-out  wings. 

The  projectila  has  canards  to  control  attitude.  Ideally  tha  projectile 
has  no  spin;  however,  due  to  effects  from  miscellaneous  foroes  and  soasnts 
which  work  on  the  body,  it  nay  acquire  spin  motion. 

In  order  to  calculate  the  trajectory  cans ids ring  such  conditions, 
solid  body  equations  with  6  degrees  of  freedom  are  induced.  One  of  the 
main  problems  which  will  occur  during  flight  motion  of  the  projectile  Is 
trajectory  errors  which  arm  caused  by  numerous  origins  especially  due 
Co  the  wind. 

Total  range  of  these  projectiles  are  rather  great.  Consequently, 
perturbing  influences  from  these  origins  arm  Important  to  notice  when 
long  range  trajectories  an  calculated.  Soma  of  tha  primary  origins  of 
error  are  mesa  unbalance  due  to  Impreciseness  of  construction  and  from 
shock  of  launching,  misalignment  of  jet,  misalignment  of  intake  air  flow, 
and  so  on. ^ 

In  order  to  minimise  these  effects,  it  is  useful  to  sat  tha  value  of 
thrust  of  the  ramjet  equal  sad  opposite  to  that  of  serodynmsic  drag. 

Tha  trajectory  is  equivalent  to  projectile  flight  in  e  vacuum  and  is 
known  as  thrust-equal-drag  trajectory.  The  effect  of  wind,  whan  tha 
direction  of  wind  is  tha  same  as  that  of  the  trajectory,  cm  be  readily 
calculated.  In  case  of  a  aide  wind,  tha  projectile  will  be  carried 
laterally. 
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ra  order  Co  can cal  the  effect  of  a  aid*  wind  which  will  carry  tha 
projectile  latarally  and  to  restore  tha  planar  notion  of  tha  trajectory, 
active  control  can  ba  used.  For  example,  canards  can  be  employed.  Due 
to  active  aerodynamic  control,  total  range  will  be  affected.  Active 
control  nodule tea  aerodynamic  drag.  The  boundary  and  initial  conditions 
will  Influence  whether  total  range  will  be  expmded  or  shortened  by  the 
effect  of  canards. 

In  this  paper,  affects  due  to  winds  are  considered;  however,  pertur¬ 
bations  of  other  origins  are  concentrated  in  one  tern  In  the  equations  of 
notion  an  constants.  In  order  to  calculate  these  data  by  electronic 
conputer,  flaw  diegrane  for  FOBSBAH  are  Indicated. 

2.  Coordinate  System 

In  order  to  construct  solid  body  equations,  coordinate  system  are 
considered  as  follows: ^ 

(1)  Ground  axle  coordinate  system  (0-xyz) . 

The  origin  la  on  the  canter  of  gravity  of  die  projectile  at  the 
dm  of  departure  of  projecdle.  The  x  axle  la  toward  the 
direction  of  horizontal  line  of  shot.  The  y  axis  is  vertically 
uprlgit,  and  the  a  axis  is  horizontally  rectangular  in  a  right 
hand  sense. 

(Z)  Velocity  sxU  coordinate  ayeten  (Q-x^y^z^) . 

Tha  origin  la  at  tha  canter  of  gravity,  CG,  of  tha  projectile. 

The  xy  axis  la  along  tha  velocity  vector  of  the  projectile,  l.a. , 
x^  is  tangent  to  the  trajectory.  The  transformation  to  create 
the  (Sysy^fSy)  coordinate  system  from  die  (x,y,s)  coordinate 
ay  a  ten  is,  first,  a  rotation  about  tha  z  axis  by  an  angle  8. 

A  new  coordinate  system  la  fozmsd  and  is  identified  as  the  (x^,y^,.z) 
coordinate  system.  . 
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Hie  velocity  vector,  v  and  axe  located  in  die  x^-z  plane. 

Nov  rotate  the  (x^,y^,z)  coordinate  systen  about  the  y^  axis 
until  the  x^  axis  is  parallel  with  v.  The  angle  of  rotation  about 
the  y^  axis  is  a.  The  resulting  coordinate  system  is  the 
(xy,  yy,  zy)  system.  Befer  to  Figure  1  for  an  illustration  of 
the  axes. 


Figure  1.  Illustration  of  geometry  for  the  velocity  axis  coordinate 
system. 
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(3)  Mon  rotating  body  axil  coordinate  system  (0  -  x^y^r^) . 

The  origin  ie  at  the  center  of  gravity  of  projectile:  the 
Xg  alia  coincides  with  the  longitudinal  axle  of  syne  try.  The 


yB  and  Zg  axis  axe  transformed  as  the  aane  wanner  as  that  of 
velocity  axis  coordinate  system  by  $  instead  of  0,  and  by  5  in¬ 
stead  of  a. 

(4)  Wind  axis  coordinate  system  (0  -  xva?va*va)»  The  system  moves  with 
relative  wind  speed  similar  to  velocity  axis  coordinate  system. 

The  relevant  angles  axe  (0^j#)  instead  of  (8,c)  with  the  x^a 
axis  along  the  relative  velocity  vector. 

1  Equations 

Solid  body  equations  with  six  degrees  of  freedom  are  developed  by 
the  method  of  coordinate  transformation  in  order  to  calculate  the  trajectory 
of  the  projectile  with  ramjet  propulsion  as  follows: 
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-1  sinO  cosC  cos(<ft-8)  -  cos a  sing 


(1>) 

(14) 

(15) 

(16) 

(17) 

(18) 

(19) 

(20) 


coso  cos£  cos(<p-9)  +  sino  sinC 
For  a  projectile  with  ramjet  propulsion,  eq.  (16)  is  true  usually.  At 
launch,  p  m  0,  and  there  is  no  initial  spin  so  that  p  »  0.  For  case  of 
p  ■  0,  equations  (11)  and  (12)  become 


q-MyB/i  (ii)’ 

r  -  MzB/I  (12)' 

Furthermore,  if  thrust  has  a  value  near  to  drag,  then  dispersion  becomes 
a  minimum  so  that  the  last  terms  of  equations  (14),  (15),  (17)  and 
(18)  are  deleted  from  these  equations. 

Effect  of  Vlad 

When  a  wind  blows,  a  wind  axis  coordinate  system  is  considered. 

Using  these  coordinates,  the  x  axis  is  along  the  relative  velocity  between 
projectile  and  air.  The  motion  of  the  projectile  can  be  calculated  the 
same  way  as  in  the  case  of  no  wind.  For  this  case,  we  use  a  subscript  (a) 
which  is  applied  to  the  notations  of  the  case  of  no  wind. 

In  order  to  obtain  velocities  which  were  observed  from  the  ground , 
we  must  add  distances  transported  by  wind  during  the  time  interval  to 
values  calculated  on  the  wind  axis  coordinate.  The  relations  between 
velocities  on  wind  axis  coordinate  and  velocity  axis  coordinate  are  as 
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(21) 

a 

a 

a  x 

V 

a 

sin  9 

a 

cos 

o  »v  sin  9  cos  o  -  w 
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V 

a 

sin  ya 
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z 

(23) 

From  Chese  aquations,  we  can  calculate  v  ,  0  ,  o  knowing  data  without 

a  S  cl 

wind,  v,  9,  a  and  the  component  of  wind  relative  to  the  ground  w  w  w  . 

x  y  z 

In  the  case  of  computer  calculation  by  a  step-by-step  method  we 
can  use  (21),  (22),  (23)  to  get  Initial  value  of  v,  6,  a  In  each  suc¬ 
cessive  step.  These  values  of  v,  9,  a  are  relative  to  the  ground  axis 
coordinates,  so  that  the  trajectory  is  obtained  from  them. 

4.  Effect  of  Canard 

From  mentioned  above.  It  is  clear  that  projectile  is  not  within  the 
plane  of  launching  but  flies  in  a  different  direction  from  the  first 
schedule  when  wind  is  blowing.  In  order  to  cancel  the  miss  distance, 
we  can  use  a  canard  giving  some  suitable  trim  angle.  Trim  angle  will  be 
given  on  body  axis  coordinate  and  whose  effect  will  act  as  sinusoidal 
form  such  as  v  sin  y  where  y  is  factor  angle  measured  from  y_  axis  right 

•D 

wise  around  x^  axis.  In  order  to  cancel  the  effect  of  sidewise  dispersion 

•  a 

z  *  0  so  a  *»  0  ,  0.  «*0  must  be  satisfied, 
therefore 

FZv  ‘  0  (24) 

In  order  to  demand  (24)  we  must  add  a  term  for  the  canard  on  the  right 
terms  of  eq.  (15) 

-mg  sin  9  sin  a  -  TS  +  *spv^SCT  6  +  C  sin  it  •  v  sin  y  (25) 

La  e 

There  are  many  factors  to  give  effect  to  the  left  side  terms  of  this 
equation,  so  it  is  difficult  to  discuss  analytically  However,  the 
value  of  the  terms  are  calculated  on  each  step  so  that  we  can  decide  a 
suitable  value  of  canard  angle  for  each  step  of  the  calculation. 
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5.  Program  Flowcharts 

In  order  to  calculate  these  equations*  a  computer  program  will 

(3N 

be  built.  It  is  illustrated  by  program  flowcharts.  This  computer  program 
consists  of  a  main  program  and  subroutines. 

Inputs  which  are  used  for  calculation  are  ready  and  stored  in  incre¬ 
ment  address  by  orders  in  the  program  to  use  for  calculation.  Data 
calculated  are  provided  printed  on  sheets  and  so  on.  Program  flowcharts 
of  main  program  and  subroutines  are  shown  on  figures*  and  they  are  explained 
as  follows; 

a.  Main  Program 

Main  program  has  a  job  to  store  and  supply  data  from  incre¬ 
ment  address  of  computer  and  to  make  Initial  condition  in  shape 
suitable  for  calculation.  It  gives  these  data  of  subroutines  to  com¬ 
pute  and  to  get  results  of  calculation.  Data  returned  from  sub¬ 
routine  are  printed  periodically.  Calculation  is  carried  by  one 
step  of  variables. 

These  data  are  accumulated  step-by-step  until  they  come 
to  some  objective  point.  When  the  objective  point  is  satisfied* 
the  calculation  is  stopped.  The  point  where  the  projectile 
touched  the  ground*  point  of  print  and  points  to  change  initial 
conditions  are  detected  and  ordered  to  computer.  It  contains  FORMAT 
to  print  these  results, 
b.  Subroutine  C01€UTE 

m  order  to  calculate  eq.  (1)  *  (12)  there  need  data  of  eq.  (13) 

*v>  (20)  which  are  in  the  state  suitable  for  input.  In  this  subroutine 
these  data  are  prepared  using  three  subroutines  of  ATOM,  FORCE  and 
MDMSM.  These  data  are  sent  to  subroutine  RKG.  There  eq.  (1)  ^  (12) 
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are  calculated  and  results  are  sent  back  to  thla  routine.  Additional 
data  are  calculated  here.  These  data  are  sent  back  to  naln  program, 
t.  Subroutine  BKG 

This  subroutine  has  a  job  to  calculate  each  of  the  dif¬ 
ferential  equations  In  a  Run ga-Kutt e-Glll  method.  Suitable  data 
for  Input  are  sent  f roe  subroutine  COMPUTE.  Thla  subroutine  calls 
differential  equations  one  by  one  fron  subroutine  DANDO. 

d.  Subroutine  DANDO 

This  subroutine  stores  each  equations  wfawch  coapoees  the  tra¬ 
jectory  equations  in  state  of  FORTRAN  fora. 

e.  ATOM 

Thla  subroutine  provides  air  teaperature  and  air  density  for 
the  standard  state  of  heights  froa  aea  level.  Air  teaperature  Is 
used  to  calculate  Mach  number. 

f.  FORCE 

Data  of  forces  are  stored  in  Increasnt  address.  On  this  sub¬ 
routine  forces  that  are  arranged  suitably  to  the  state  of  input  of 
the  instance  are  selected  and  prepared  for  supply.  This  subroutine 
has  a  function  to  asks  linear  interpolation. 

g.  MOWN 

Data  of  momenta  are  stored  in  increasnt  address.  Moments  that 
are  arranged  suitably  to  the  state  of  Input  of  the  instance  are 
selected  and  prepared  to  simply.  This  subroutine  has  a  function 
to  make  Unaar  interpolation. 
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Progr 


Flow  Charts 


(a)  Main  Program 


Print 


6.  Conclusions 


Solid  body  squstlons  with  6-degree-of-f resdosi  are  formulated  in  order 
to  calculate  the  trajectory  of  a  projectile  which  la  flying  vith  the  thrust 
of  a  ramjet. 

In  order  to  calculate  by  electronic  digital  conputer;  program  flow 
charts  are  developed.  FOftTBAH  programs  will  be  made  for  each  subroutine  and 
program.  Trajectory  errors  due  to  die  effect  of  the  wind  are  considered 
mainly.  Dispersion  due  to  the  wind  will  be  calculated  by  these  program. 

If  there  need  to  cancel  the  effect  of  a  side  wind,  trim  angle  of  canard  will 
be  given  step  by  step  in  the  procedure  of  calculation. 
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